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Abstract. An inverse polynomial has a Chebyshev series expansion 

k oo 

V^bjTjix) = ^'a„T„(x) 

j=0 n-0 

if the polynomial has no roots in [—1, 1]. If the inverse polynomial is decom- 
posed into partial fractions, the a n are linear combinations of simple functions 
of the polynomial roots. Also, if the first k of the coefficients a n are known, the 
others become linear combinations of these with expansion coefficients derived 
recursively from the bj's. On a closely related theme, finding a polynomial 
with minimum relative error towards a given f(x) is approximately equiva- 
lent to finding the bj in f(x)/ J^* bjTj(x) = 1 + X]fe+i a n T n {x), and may be 
handled with a Newton method providing the Chebyshev expansion of f(x) is 
known. 



1. Introduction and Scope 

The Chebyshev polynomials T n (x) are even or odd functions of x defined as |2 
(22.3.6)]0 (3-6)] 

(1.1) T Q ( X ) = 1, T^^i-r^^i^)^, , = 1,2,3... 

m=0 v ; 

where the Gauss bracket |_-J denotes the largest integer not greater than the number 
it embraces. The reverse formula is X, p. 412] |13l p. 52] |14| 

where the prime at the sum symbol means the first term (at j = and even n) is 
to be halved. The polynomials are orthogonal over the interval [—1, 1] with weight 
function l/\/I - x 2 [TJ (22.2.4)] 7, (4.2)] 



dx 



7r , n = m = 0, 



(1.3) / T n {x)T m {x)-== = <^ vr/2, n = m ? 0, 



vT 



, n ^ m. 
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The product rule is [i, (22.7.24)] [El (A.1)]HH (2.10)] 

(1.4) T n (x)T m (x) = - (T\ m _ n \{x) + T m+n (x)) . 

The indefinite integral is (4.8)] [13 p. 54] [U (2.12)] 

Ti{x), n=0, 

n = 1, 

n > 1, 



(1.5) / T n (x)dx = { 

' 1 T n+1 {x) _ T„-i(a) 
2 ^ n+1 n-l 



2 

which correlates to the derivative 



d 



(1.6) —T n {x) = 2n J2 T '( 



dx 

1 = 
n-l odd 



The expansion of an inverse polynomial of degree k in a power series is 
(3.6.16)] 



(L7) r^7 = £ 

with recursively accessible 0.313] 

n-l 

-j<k 



1 1 

(1.8) c n = -— y~] dn-jc-j, c = — , n > 1. 

"o ~£ d 



The topic of this script is the equivalent arithmetic expansion of the inverse 
polynomial in a Chebyshev series, 

oo 

1 1 v-V 

Ej 'I J- 1 '-' L,j=O b 3 T 3( X ) n=0 

i.e., computation of the coefficients 

ft irA 2 f 1 Tn ^ dx 

(1.10) a, — ' 



given the sets {bj} or {dj} that define the original function. Both sets are related 
via [S3 (3)] [SI (37)] and with JTj) via 

oZ-l k (2+1 — 1 V 

(1-n) * = V S (_) iL fesr^ ' = °-- fe - 

l 2 /' 



The expansion l|1.9|) exists if the inverse polynomial is bound in the interval 
[—1,1], i.e., if ~^2,djX 3 has no roots in [—1,1]. 

Characteristic generic methods of evaluating (|1.10|) are not reviewed here: (i) 
Fourier transform methods [7J f4.7')][5 HT2*lllO) . (ii) sampling with Gauss-type quadra- 
tures ( 25.4.38')] |25U33ll2U) which effectively means using an implicit intermediate 
interpolatory polynomial to represent 1 / X^=o ^j'^j ( x )i ( iii) approximation by trun- 
cation of l|1.7f) . then insertion of (|1.2|l . (iv) using the near-minimax properties of 
the Chebyshev series [251 12T] . 
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Remark 1.1. The Fourier-Chebyshev series [EI EZ| [2HI (3.4.1f)] 



(1.12) Tn{z) tT \;-™\^ = y Tkm+n {z)t k 

V ' l-2tT m (z)+t 2 ^ + V ' 

my 1 fc=0 

provides special cases of polynomials with particularly simple expansions. 
Remark 1.2. Insertion of n = 1 in l|1.4l) shows that the coefficients of 

oo 

(1.13) f(x)=J2'fMx) 
and 



(1-14) lM = J2'g n T n (x) 

n=0 



are related as 

(1-15) fa = gi, 2/ n _i = g n _ 2 + g n , n > 2. 

Chapter[21explains how the a n of (|1.1U|) could be computed supposed the inverse 
polynomial has been decomposed into partial fractions. Chapter[3]provides a recur- 
sive algorithm to derive high-indexed a n (n > k) supposed the low-indexed a n are 
given by other means. Chapter 0] touches on a (standard) integral- free method to 
compute approximate low- indexed a n , and Chapter deals with a specific inverse 
problem — which is finding the bj from partially known a n — related to polynomial 
approximants with minimum relative error. 

2. The Case of Known Partial Fractions 

The straight way of computing the Chebyshev series uses the decomposition of 
lfYldjXp into partial fractions 2.102], which reduces (|1.9I) to the calculation 
of the a n ^ s in 

(2.1) 1 = a n , s (z)T n (x), 

^ ' n=0 

where z is a root of the polynomial. 

fc 

(2.2) X''-'-' "• 

3=0 

Sign flips of z and x in 1|2.1|) show that 

(2.3) a n , s (-z) = (-) n+s a n Az). 

The case of s = 1 has been evaluated earlier ^| (A. 6)] [23 02 based on ^ 
(22.9.9)]|33 (18)], 

(2.4) a nA {z)= — w = z+(z 2 -l) 1 ' 2 , zi[-l,l\. 
The branch cuts of (z 2 — l) 1 / 2 must be chosen such that > 1. 
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Example 2.1. 

(25) = 

y ' ' l + x 2 2i-x 2-i-x 

consists of two terms, 

(2-6) a n i(i) = — — , a n i(-i) = (-) n = — , 

V ' ' U ' (l + %/2)« ' V ' V ' {l + s/2) n 

which recombine with the two factors i/2 and —i/2 to [2S1 (3.4.1a)] 

(2.7) s = \/2 V 1 ; _ T n {x). 

l + x 2 ^ (1 + V2) n 

71=0,2,4,6,... v ' v ' 

Remark 2.2. The shifted Chebyshev polynomials T*(x) = T(2x — 1) are orthogonal 
over [0, 1] with weight l/^x(l - x) 1, (22.2.8)] [2l>|. From we get 

OO 

(2.8) ^_^ = 2 ^' a „, s (2z-l)T„*(x), 
and l|1.5[) becomes 

^T(x), n = 0, 

(2.9) /T:(x)d,= <| ,^ >w t ,_ iWN »=L 



Example 2.3. An example of s = 1, z = — 1 in l|2.8f) is 

OO 

(2-10) -J- = -— — = -2j^ a n<1 (Smx), 

l + x — 1 — x 

n=0 

where 

( _ r+1 



(2.11) a „i(-3)= _ 

1 ' ' U ; x/2(3 + 2V2)» 

according to l|2.4|l . 

Higher second indices s of the a n , s are obtained from l|2.1|l by repeated derivation 
w.r.t. z, 

via [13 0.432.1], 

r9i^ o L y (-)' A A (2^)- 2 ' >n 

with Pochhammer's Symbol defined as (6.1.22)] 

(2.14) (a) k = a(a + l)(a + 2)---(a + k-l)=T(a + k)/T(a), (a) = 1. 
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The formula 

(S - 1) / - — = r r + A s 

K 'J {z-x) s (z-xy- 1 



GO GO 

(s - 1)^ a„ iS / T n (x)dx = ^2 a n>s -iT n (x) + A s , s>2 



n=Q " n=0 

in conjunction with the method quoted by Cody Q (4.8)] [23 (25)] yields 

2n 

(2.15) a n+ \ Jz) = a„-iJz) -a n s -x(z), n > 1, s > 2. 

s — 1 

One needs l|2.13() and 

2 f 1 T x (x) dx 

(2.16) = -ao, s (z) + za , s+ i(2;) 

to start the recurrence 12.15|l and to obtain all coefficients in 1|2.1[1 for a particu- 
lar z. Closed form expressions for solving these recurrences in terms of Legendre 
Polynomials of z/y/z 2 ~ 1 = (w 2 + \)/(w 2 — 1) have been given by Elliott [B]. 

Remark 2.4. (|2.16|l may be generalized to 

(2-17) i t Xl dx = "t(-r(!)z l - m ao,n-m, Kn. 

ir J_ x {z - x) n ^J\-x 2 ^ \mj 

and with (|1.2|l and (|1.4|) to 

r~ o1 2 f 1 x l TJx) , 1 J^> ( l\ r n 

(2 - 18) vi -^3 h yh- S i,+^]- 

Example 2.5. An example of degree k = 3 is 

1 1 

(2-19) 



(4 - x) 2 (5 + x) 78±Tb(z) - 23±Ti(a;) - l\T 2 {x) + ±T 3 (:r) 

111111 
( ' ' ~ 9 (4 - .t) 2 + 81 (4 - x) ~ 81 (-5 - x) ' 

The root at z = 4 yields 

(2.21) a ,i(4) = 2/V15 w 0.5164 

from and 

1 2-4 

from (|2.13|) . The root at z = — 5 yields 
(2.23) a ,i(-5) = 2/(-v / 24) » -0.4082 

from (|2.4(l . The combined total in (|2.20j) is 

(2.24) 



(2.22) ao2 ( 4 ) = 2 .-.-— 3 « 0.1377 



a = - [ - -i- - dx w --0.1377+— -0.5164- — -(-0.4082) w 0.0267. 

7r (4 — x) (5 + x) \f\^x 2 9 81 81 v 
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Example 2.6. A case of k = oo is P3 1.421.2] 



(2.25) 



tanh(7rx/2) 



1 



x tt^2(2to-1) i(2m-l) 

The roots at z — ±i(2m — 1) yield 

(2.26) (In. 



-i(2m — 1) — x 



A (z) = 2/ (±i-\/4m 2 -4m + 2) , 



and the combined total is 
(2.27) 

2 



a 



1 tanh(-7ra;/2) dx 
-1 x 



E 



Vi - x 2 7T ^ (2m - ^^m 2 -m + 1/2 



2.38. 



Remark 2.7. From l|2.4() 

9a„ 1 (z) 

(2.28) 



_Z 2 -1 (z 2 -l)!/2_ 

so the (linear) propagation of the absolute relative error in the root z to the error 
in the coefficient a rii i is 

(2.29) 



Aa„ 1 



1 



(z 2 - l) 1 / 2 



Remark 2.8. An associated factorization 53j=o bjTj(x) cx IlL^iX 2 ™ — x) Sm , with Z 
different roots of multiplicities s m , decomposes the square root of the polynomial 
into a Z-fold product of series of the prototypical forms 



(2.30) 



z — x 



^ q n (z)T n (x), s m = 1, 
zT (x) - Ti(x), s m = 2, 



(2.31) 
where Q2I 2.576.2] 

2 r* 4 

(2.32) 170(2) = - / dWz-cost = ~\Jl + zE{- 

V Jo 7T 1 



+ z 



is related to Complete Elliptic Integrals of the Second Kind E in the notation of 
(17.3.4)]. The q n (z) with n > 1 follow recursively using [TJ (17.1.4)]. In particular, 
one may expand T n {x) in terms of P^' 1 ' 2 \x) with ^2 (1-4)] to obtain 



(2.33) 



?n(l) 



2 5/2 y> (-n)i(n)j 



t — (3/2),(l/2)r 



?i = 0,1,2, 



for the Chebyshev coefficients of \J\ — x. See for an application. 

3. Recurrence of Expansion Coefficients 

The T n in QLlOj l may be decomposed into a unique product of a polynomial by 
the denominator plus a remainder of polynomial degree less than k. [The depen- 
dence on x is omitted at all T„(x) for brevity.] 



T n = (d n) T + drT 1 + --- + d^ k T n - k )(b T + b 1 T 1 + --- + b k T k ) 

(«) 

H — ;p J o + Ci Ji + c 2 J 2 + • • • + c k _ 1 lk-i- 



Mr 



An) 



(3.1) 
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Expansion with (|1.4J) yields a system of linear equations for the vector of the un- 
know 
(3.2) 



knowns dj and cj n 



V o 



/ 1 






1 


. 


. . 0| 
. . 


2b 
bi 


bi 

b + b f 


61+63 62+64 
2 2 











l>2 


bi +b 3 
2 


U 1 64 61+65 
b + 2 2 



1 



... 


... l|6 fc _i 


bk-2+bk 
2 


bk-3 
2 


... 


I) b k 


b k -l 
2 


bk-2 
2 




:| 


bk 
2 




bk-i 
2 



"fc-3 

2 



bu- 



\ 



bk 

2 / 



( ^ \ 



J") 
c fc-l 

"0 

d (n) ' 

"n-fc-1 

v 4 n i y 



The (n + 1) x (n + 1) cocmcicnt matrix A r ^ c (row index r and column index c from 
to n) is an upper triangular matrix. It hosts a k x fc unit matrix in the upper 
left corner, and is symmetric w.r.t. the minor diagonal that stretches from Aq^ to 

Aji—k.n • 



(3.3) 



A r ,c = $r,c, < C < k — 1. 



(3.4) 



.4 



.4 



c,fc+r 




r = c = 
r = 0, 1 < c < fc 

r = c= 1, 2, . . . , n — k 



This works with the auxiliary definition 

(3.5) bi = 0, i > fc or i < 0. 

Insertion of (|3.1|) into (|1.10() yields 



fc-i 



(3.6) 



2dW+ y'f^Oi, „ /.-. 



i=0 



which means that entire sequence a„ can be generated recursively from its first k 
terms, if the d Q n ^ and & are generated at the same time via (I3.2H or an equivalent 
method. Iterated full solution of l|3.2l) can be avoided through recursive generation 
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of the set K (n+1) ,c 4 ( " +1) } from {dj n \e£ n) } and {^ n_1) ,c i (n - 1) } as follows: 

(3.7) 
(3.8) 

(3.9) 4" +1) 

c 



(3.10) 



2 

>+i) 





c fc-l 
h 


d (n-l) 






2d { n) - 


f 4" } - 


d (n-l) 

"1 ! 








a j + l 




j = 2,3, . fc- 


- 1. 


c i 


u (™) 








bk 


2 : 






(n) 
C i-1 "> 


(n) 
' C i+1 " 


b k 


c?-^, J = 1,2,.. 





(3.11) Cj 
where the auxiliary definitions 

(3.12) 4" )=0 ' or J <0 ' 

(3.13) = °' j>n-k, or j < 0, 
are made to condense the notation. 

Proof. Multiply (jSTTJ by 2Ti and use (JOJl as 

n— fe 

2T 1 ^4" )r j = 4" )T o + (24 n) +4 n) ) T i 

n—k— 1 

(3-14) + £ (d™ + 4" ) 1 )T 3 + d^T*.* + ^fcTn-k+i, 

4=2 

(3.15) 2T! c^T, = c^T Q + + cf^T, + c k % T fe _ a + ^T, 

4=0 4=1 
Rewrite the last term in the previous equation 

(n) fe (n) (n) 

(3.16) C jW T* = ^IV^ - %i6oT - ... - ^i^.^x. 

&fe frT h b k 

Construct 

r («) 

2TiT„ = [(4") + ^-i)r + (24" ) +4" ) )T 1 

4=2 j=0 

(n) fc-2 (n) (n) 

+(c\ - -g—b )T Q + 2^(S-i + c }+i - -^— b i) T 3 + ( c k-2 - -j—b k -i)T k -i, 

k . =l k k 

and subtract T„_i for identification of the 4" an d c j" +1 ''> 

n-fc+1 k 

(3.17) T n+1 = 2T 1 T, i -T„_ 1 = ( ^ 4" +1) T,-)(^6 i T i )+^'4 n+1) T :; , 

4=0 4=0 3=0 
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□ 



Example 3.1. For l|2.19(l . we obviously have 



(3.18) 



.(1) _ „(1) _ „(2) _ J2) _ 



„\ L > — — n \") — n J. 1 ) — r-^ — 1 

L 2 ~ L ~ L l ~ u i L l — L 2 ~ 1 



in lpTT]l , The formulas J3~7| - ((3~TT]| predict at n = 2 
(3.19) 



7 (3) 



1 



,(3) 



1% \ ' 1 J3) 



1 



-23i • 1 



1, 



,(3) 



1/4 



1/4' 2 1/4 ' ~ x ' 1/4 

With these, (|3.6j) gives at n = 3 

(3.20) a 3 = 8+(-314)-a + 93-ai+6-a 2 
which is correct since 

(3.21) a w 0.02671606, m « 0.00412578, a 2 w 0.00087916, a 3 w 0.00013030. 

The next step of the recursion is 
(3.22) 



6 



cn 4) 78^-6 (4) , -23]-6 

^ = 93-^, ^)= 2 .(-314)+6-^^ > 



C W - 93 12— 

C 2 1/4 



-I7 



(3.23) a 4 = 48 + (-1791) • a + (-64) • a x + 128 • a 2 
which is also correct with 

(3.24) a 4 « 0.00002159. 

4. Approximation by the Truncated Chebyshev Series 

Approximations a n to the a n of l|1.9|) may be calculated assuming that the a n 
are negligible beyond some index TV: 

AT 

(4.1) 1 n£a n T n (x). 

Y,j=0 b 3 T i( X > n=0 

If this equation is multiplied by 2 bjTj , and we stay with l|3.5(l to keep the notation 
simple, 

(4.2) 



N i k+n 

2 ~ °" n 22( b n-l + b l+n + k-n)Tl. 

n=0 1=0 



If the coefficients in front of To to Tjv are set equal on both sides, 
equations for the a n ensues: 



/ b 



(4.3) 



&2 h 

bi 2b + b 2 bi + b 3 b 2 + b± 

&2 bi + b 3 2b + 64 61 + b 5 

b 3 h + &4 h + b 5 2b + b 6 



\ 



( a \ 



a system of linear 
/ 2 \ 



/ \ a N J 



V J 
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(4.4) B r . c = 



where the coefficient matrix B r>c is symmetric and has a band width of 2k + 1: 

b c , r = 

b r , c = 

2b + b 2r , r = c^0 

6| r _ c |+6 r+c , r^c, c > 0, r>0 

This gives access to a set of approximate, low-indexed a n with no need to evaluate 
integrals nor reference to the roots of bjTj. 

Example 4.1. Again for 1|2.19[) . the choice of N = 3 yields the coefficient vector 

(4.5) a = 0.02671602, a ± = 0.00412567, a 2 = 0.00087845, a 3 = 0.00012696. 

At N = 4, this improves to 

a = 0.02671606, a x = 0.00412578, a 2 = 0.00087914, a 3 = 0.00013019, 
^■ b > a 4 = 0.00002111, 

which is close to the exact results in (I3.21|) and 13.2411 . At N = 5, this improves 
further to 

, r a = 0.02671606, a x = 0.00412578, a 2 = 0.00087916, a 3 = 0.00013029, 
^■') a 4 = 0.00002158, . . . 

Remark 4.2. The matrix in (|4.3|l is the approximate, square upper left (N + 1) x 
(N+l) submatrix of the "exact" solution. The approximate solution obtained could 
be considered as if the terms Y^n=N+i a n{b n -i+b n +i + bi- n ) of I4.2fl in the i'th row 
of the system of linear equations had been neglected (as if the columns N + 1 up to 
I + k had been chopped off) . The neglected sum is nonzero only if / > N + 1 — k. An 
idea of an improvement of this algorithm is: reduce the a n (N+l < n < l + k) in the 
neglected terms via (|3.6ll to a linear combination of ai,...,fe, and re-introduce (add) 
these components (in)to the matrix — add the constant to the right hand side — in 
these rows / > N + 1 — k. This update of the system of linear equations reduces the 
rank of the matrix and does therefore not improve on what is obtained from l|4.3|l . 

The algorithm may be extended to the division problem of finding the a n from 
given /„ in 

f(x) N > 

(4-8) — k ■ ' w ]T a n T n (x); f(x) = £ f n T n (x), 

l^j=objlj( x ) n =0 n=0 

with the right hand side in l|4.3(l replaced as follows: 



(4.9) 



' fr, r = 0, 
2f r , r= 1,2,3, 



B TtC a c — < 

r— n ^ 



Example 4.3. The Chebyshev series of f(x) = sin(^x)/x starts with [5ll3l7| 



! l ' 7T ^ 2S ~, 2.552557924804531760415274, 



<«»» A = 'gp^iW(j) 

(4.11) / 2 w -0.2852615691810360095702941, 

(4.12) U ~ 0.009118016006651802497767923, 

(4.13) f e w -0.0001365875135419666724364765, 

(4.14) / 8 w 0.000001184961857661690108290062 
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(415) f = M(-) n/2 Er=n/ 2 4 + iW2), n even, 

\ , n odd. 

(Schonfelder lists 2/2 n /7r for n < 16.) If we approximate f(x) by the polynomial 

.4 

X) n =o fnT n (x), calculation of the a ra in 

fC 1 N 
(4.16) /W « V a„T n (x) 



F,' i=o fi T i( x ) 



via (|4.9|l at A~ = 8 predicts the relative error 

iV 

^'a„T„(a;)-l » -6.74 • l(r 8 T (a;) - 9.97 • 10" 7 T 2 (a;) 

n=0 

(4.17) -1.23 • lO^T^rr) - 1.09 ■ 10" 4 T 6 (a;) - 1.13 ■ lCT 5 T 8 (x). 

Remark 4.4. The functional relations l|4.8|) hold also for the shifted Chebyshev 
polynomials T* (x) : 

(4 - 18) ^ TLt ^ J2'^T:( X y, m^Y.'f^ix), 

l^j^Vj 1 j ( x ) n=0 n=0 

5. Chebyshev Approximation for the Relative Error 

The previous example of a truncated Chebyshev series had a maximum absolute 
error estimated at Y^n=6 ~ 0.000138 if terms up to k — 4 were retained, and the 
maximum relative error of the same polynomial was estimated at 53'n=o l^n| ^ 1 ~ 
0.000134 — dominated by the &e term in 14.17JI . To optimize the approximation 
of f(x) for the relative error in [—1,1], one would rather like to find the k + 1 
coefficients bj in l|4.8|l which force the relative error to be close to zero in the sense 
of 

(5.1) a — 2, a% = a 2 = a 3 = . . . = a k = 0. 

As an inversion of the problem of Sec.^J the matrix B in l|4.9|l is presumed unknown 
(up to some symmetry), and the first fc+1 elements of the vector a c and all elements 
of f r are known. The rationale is that removal of the ripples of T\{x) to Tk(x) from 
the quotient expansion leaves a quotient with an appropriate number of "critical" 
points required by the alternating maximum theorem (3EHSE1- 

Remark 5.1. The case r = in (|4.9|l in conjunction with l|5.1|l mandate 

(5.2) b = fo/2. 

Finding the constituents bj of B that solve the bi-linear equation (|4.9|) may 
proceed with a vectorized first-order Newton method as follows: 

• Chose a start solution bj, for example the obvious 

' /o/2, j = 
fj, j = l,2,. ..,k 

• Compute the a n (n — 0, . . . , A^) from bj by solving the linear system of 
equations (|4.9(l . 



(5.3) b 
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(5.4) 



(5.5) 



Compute an approximate (N + 1) X k Jacobi matrix 



J r.c 



by partial derivation of the first N + 1 equations of l|4.9|l w.r.t. the bj, i.e., 
by solving the k systems of N + 1 linear equations 







dao 


da 




,)l n 




■ db k 




Da j 


aa\ 


oai 






9&i 


db 2 


■ db k 




V 


daN 


daN 




) 


dbi 


db 2 


■ db k 



( 



N 



«i 



a 2 



«3 



a + a 2 Oi + a 3 a 2 + 04 
d\ + a 3 a + 64 d\ + a 5 



flJV-2 



dN-3 



flfc-1 



OAT-fe-l 



flfc-1 + flfe+1 



O-N-k 



J 



for r = 0, . . . , N and j = 0, . . . , k — 1. The column da c /dbo of the Jacobi 
matrix is not calculated, as 60 is assumed fixed according to l|5.2|) . 
• Compute the next iterated solution &j + Aj (j — 1,2, ... , fc) of the polyno- 
mial coefficients by solving the system of k linear equations 



(5.6) 



^ db, 3 



3 = 1 



-at, 



1 = 1, 



for the first-order differences Aj. This equation is the first-order Taylor 
expansion of di as a function of the bj set to the target (|5.1|) for this update. 
The k x k coefficient matrix dai/dbj is a square submatrix of the Jacobi 
matrix calculated in the previous step. 
• Return to the second bullet for the next cycle until the do to dk are suffi- 
ciently close to l|5.1|). 

Remark 5.2. This algorithm involves only /o to /jv, but no higher order approx- 
imants to f(x). It therefore adapts a polynomial of degree k to a polynomial of 
degree N. 

Example 5.3. The error terms (|4.17|l for the polynomial Ylj=o bjTj (x) change to 

N 

J2'a n T n {x)-l pa 5.2- 1(T 12 To(.t) +4.7- 10" n r 2 (x) +6.3- 10" 12 T 4 (a;) 



(5.7) 



-1.08 • 1(T 4 T 6 (2;) - 1.11 ■ 10- 5 T 8 (x) 



after one Newton iteration, reducing the relative error to J2'n=o l^nj — 1 pa 0.000119. 
During further iteration cycles the relative error stays about the same because it is 
dominated by d^T^x) which is out of reach of the polynomial base with k = 4. 

Example 5.4. An IEEE "single" precision accuracy of f(x) = sm(^x)/x with 
a relative error smaller than 2~ 24 ps 6.0 • 10~ 8 needs k = 8. Truncation of the 
Chebyshev series for f(x) after k — 8 yields an estimated maximum absolute error 



of |/n| 

polynomial is also 



6.7 ■ 10~ 9 evaluated at N 

N , . 



n=0 



On 



1 w 6.7- 



16. The relative error of the same 
10~ 9 . After four Newton iterations, this 
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value drops to 5.9 • 10 9 with coefficients given in the following table — remaining 
very close to those cited after (|4.10l : 

_n K 

1.276278962402265880207637 

2 -0.2852615691810328617761446 

4 0.9118016006289075331306166-10- 2 

6 -0.1365874893444115901818408-10- 3 

8 0. 1184206224108742454613850- 10 -5 

Example 5.5. g(x) — cos(^x) has the expansion coefficients [51 1231 l3T)| 

2(-) n/2 J„(7r/2), neven, 
, n odd. 



(5.8) g n 



The approximation g(x) « Y^'n=o 9nT n {x) has an estimated maximum absolute 
error of Yln=k+l ~ ^-7 ' ^ 8 ^ or tne polynomial of degree fc = 8 evaluated at 
N = 16. Because g(x) is zero at both ends of the interval [—1,1], the algorithm 
does not find polynomials J2j=o bjTj(x) with a uniformly convergent Chebyshev 
expansion of the relative error — any a n obtained depend strongly on N. We there- 
fore "lift" both zeros by looking at f(x) — cos(^x)/(l — a; 2 ) instead, which has the 
expansion coefficients 

(5.9) f a = 7rJi(7r/2), 

(5.10) h = /o-2$o, 

(5.11) /„ = 2/„_ 2 - fn-4 - 4.g„_2, n = 4, 6, 8, . . . , 

(5.12) /„ = 0, nodd. 

Truncation of the Chebyshev series for f(x) after k = 4 yields an estimated maxi- 
mum absolute error of E„ =fe +i |/n| ~ 2.7 • 10" 5 evaluated at N = 8. The relative 

error of the same polynomial is Y^'n=o l°«J — 1 ~ 3.3 ■ 10~ 5 . After four Newton 
iterations, this value drops to 3.1 • 10~ 5 with coefficients given in the following table: 
_n bn 

0.8903651967922106931461297 

2 -0.1072744347398521266520654 

4 0.002332103968386755210894198 

Example 5.6. The coefficients of the Chebyshev series of arcsinoi and (arcsin x)/a; 
(App. are slowly descending. The infinite slope of arcsinoi at x — ±1 ren- 
ders both series inefficient, so we turn to iarcsin-^ instead as configured in 
HP. 10(1 . Keeping terms up to J^q yields an estimated maximum absolute error 
of E„ =38 \fn I ~ 1-4 • 10~ 17 evaluated at N = 108. The relative error of the same 
polynomial is w 1.9 • 10~ 17 . After four Newton iterations, this value drops only 
slightly to 1.8 • 10 -17 ; obviously, there is not much room to improve the polynomial 
representation w.r.t. an optimized relative error in cases where the amplitude of the 
function is small over the x-interval. 

Example 5.7. The expansion for exp(ir) in — 1 < x < 1 reads ^ (9.6.19)] [2E1 
(3.4.1e)]QJ p. 69] (33)] 

OO 

(5.13) exp(x) = 2 ^' /„(l)T„(x). 

71=0 
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Truncation after k = 14 yields an estimated maximum absolute error of X)^=fe+i l/nl 
4.9 • 1CP 17 evaluated at N = 42. The relative error of the same polynomial is 
12' n =0 I'M — 1 ~ 1-3 • 10~ 16 . After four Newton iterations, this value drops to 
7.5 • lO" 17 



T N d x 1 



with coefficients given in the following table, also listed as f{x) 
K d n 




1 

2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 



1.2660658777520083355982446 
1.1303182079849700544153921 
0.2714953395340765623657051 
0.4433684984866380495257150-10- 1 
0.5474240442093732650276168-10~ 2 
0.5429263119139437503621352-10- 3 
0.4497732295429514665443872-10~ 4 
0.319843646240199050133412M0- 5 
0. 1992124806672795001043316- 10~ 6 
0.1103677172551632915777862-10- 7 
0.5505896079551881657982078-10~ 9 
0.2497956604792065959497342-10- 10 
0.1039151254481832513826561-10- 11 
0.3990676874210170341122722-10- 13 
0.1400237499722866786358850-10- 14 



1.00000000000000002107745526254 
1 .00000000000000063548946139343 
0.499999999999997953936666685291 
0.1666666666666422610320391 

-l 



0.4166666666669875817272051 
0.8333333333602639662588442 
0. 1388888888702869286166025 
0.1984126971086418099245159 
0.2480158780231612103680909 
0.2755735152373104259316644 
0.2755725369287090362239172 
0.2504783672757589754944252 
0.2088034159586738951818317 
0.1634581247676485771723867 
0.1147074559772972471385170-10- 10 



10 

10- 2 

io- 2 

lO" 3 
10 
10 
10 
10 
10~ 8 
10" 9 



Example 5.8. The expansion for J (|x) in — 1 < x < 1 reads (El 6.681.5] 



(5.14) 



M-x) = 2 ^ 
n=0,2,4, 



(-T /2 J 2 n/ 2 (7)Tn(x). 



•4- 

Truncation after Tiq(x) yields an estimated maximum absolute error of Xm=fe+i IM ' 
7.3 • 10 -19 evaluated at N = 48. The relative error of the same polynomial is 

S ra =o I'M ~ 1 ~ 1-6 • 10~ 18 . After four Newton iterations, this value drops to 
10 -18 with coefficients given in the following table: 

n b n d„ 



1.3 ■ 





2 

4 

6 

8 

10 

12 

14 

16 



0.7252769164405135618043045 
-0.2638108118461404734713153 
0.1072184541022420669256084-10- 1 
-0.188568764213595296719917M0- 3 
0.1845983728936489887451460-10- 5 
-0.1150537142155094251800350-10- 7 
0.4965029850154789447530764-10- 10 
-0.1571252252452718608949964-10- 11 
0.3800986508122698831881511-10- 15 



A set of bj in 
(5.15) 



R(x) 



fix) 



E,=0 b T ji x ) 



0.9999999999999999991311745 
-0.6168502750680847778603892 
0.9512606546288948620024320-10- 1 
-0.6519837738512518004083602-10- 2 
0.2513602312234872245916252-10~ 3 
-0.6202064609606906421245435-10- 5 
0.1062698637612363296679714-10~ 6 
-0.1336990135568532922581048-10- 8 
0.1245507258981645953230933-10- 10 



1 



found that way is also a starting point to calculate the solution with the minimax 
property of the relative error: This locates the local minima and maxima of R(x), 
computes the mean of their absolute values, and iteratively adjusts the bj such 
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that the absolute values of the new alternating extrema equal that mean. The 
corrections Aj to the bj can be computed by expansion of (|5.15|) to first order in 
Aj keeping the abscissa of the extrema fixed, which ends up in a linear system of 
equations for the A,,-. 

Example 5.9. An IEEE "double" precision accuracy of f(x) = sin(-|.T)/a; with 
a relative error smaller than 2" 53 w 1.1 • 10 -16 needs k — 16. Truncation of 
the Chebyshev series of Example l|4.3[l for f(x) after k = 16 yields an estimated 
maximum absolute error of En=fc+i \fn\ ~ 4.1 • 1(T 19 evaluated at N = 32. The 

relative error of the same polynomial is X)' n =o — 1 ~ 3-8 • 1CP 19 . After four 
Newton iterations, this value drops to 3.5 ■ 10~ 19 with coefficients b n given in the 
following table: 



n 


K 











1.276278962402265880207637 




1.5707963267948966188688195 




2 


-0.2852615691810360095702941 




-0.6459640975062461962319336 




4 


0.9118016006651802497767923-10- 


-2 


0.7969262624616554097627533- 10" 


-l 


G 


-0.1365875135419666724364765-10 


-3 


-0.4681754135303468240882506-10 


-2 


8 


0.1184961857661690108288872-10- 


-5 


0.1604411847100114088031881-10- 


-3 


10 


-0.6702791603827441081706121-10 


-8 


-0.3598843013917326159520456-10 


-5 


12 


0.2667278599017903283863443-10- 


10 


0.5692135656122429901944357-10- 


-7 


14 


-0.7872922004615709018594325-10- 


-13 


-0.6684369436484103757933363-10 


-9 


16 


0.1791929094718284072119916-10" 


15 


0.5871793257572873247522307-10" 


11 



The actual relative error of this approximation is shown in Fig. ^ as a continuous 
line, with a maximum of 2.9- 10- 19 . The dashed line with a relative error of 2. 6-10- 19 
in comparison results from further minimax optimization with coefficients shown 
in the next table: 



3 I r 




-3 1 1 1 1 1 1 1 1 1 1 1 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

x 



FIGURE 1. The relative error R(x) to f(x) = sm(nx/2)/x for both 
tabulated parameter sets of X^nLo 2 b n T n (x) of Example 15.91 
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_n K dn 

1.2762789624022658802075437 1.5707963267948966188314659 

2 -0.2852615691810360095705230 - 0.6459640975062461915471363 

4 0.9118016006651802497528156-10- 2 0.7969262624616544421893744-10- 1 

6 -0.1365875135419666726405733-10" 3 -0.4681754135302704719117724-10- 2 

8 0.1184961857661689920542732-10~ 5 0.1604411847070460989830944-10~ 3 

10 -0.6702791603827612608171959-10" 8 -0.3598843007652658555526627-10- 5 

12 0.2667278599019855592489579-10~ 10 0.5692134921914455833455723-10" 7 

14 -0.7872921659616258733890169-10" 13 -0.6684324580312975131354658 TO" 9 

16 0.1791589025538146793760922-10- 15 0.5870678918883399413795788-1Q- 11 



Example 5.10. As an example for l|4.18|l . consider exp(a;) = Y^'n°=o f n ^n( x ) over 
< x < 1 (4.2.48)] C2]. The /„ are represented via Q] (9.6.26] through 
modified Bessel Functions /„, 



(5.16) f n = 2^/„(l/2); f n+1 = -Anf n + /„_!, 



ft I 



3.506775308754180791443893 

1 0.8503916537808109665352350 

2 0.1052086936309369253029528 

3 0.008722104733315564111612874 

4 0.0005434368311501559635982758 

5 0.00002711543491306869404046064 

Truncation of the Chebyshev series for f(x) after k = 3 yields an estimated max- 
imum absolute error 

of E„U+i \fn\ » 5.7 • 10~ 4 evaluated at N = 9. The relative 

error of the same polynomial is 2J „ = o \®n\ — 1 ~ 5.1 • 10 . After four Newton 
iterations, this value drops to 4.0 T0~ 4 with coefficients given in the following table: 
_n K 

1.753387654377090395721946 

1 0.8503902561425088936327743 

2 0.1051918520893768747555014 

3 0.008587089960927766771654559 

If we proceed to k — 12 at N = 36, the estimated maximum relative error 
becomes 6.1 • 10 -18 with the following coefficients: 
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n 



1.7533876543770903957219464 

1 0.8503916537808109665352350 

2 0.1052086936309369253029528 

3 0.8722104733315564111612874-10- 

4 0.5434368311501559635982758-10- 

5 0.2711543491306869404045765-10- 

6 0.1128132888782082788967416-10- 

7 0.4024558229870710027066467- 10" 

8 0.1256584418283842256517024-10- 

9 0.3488091362080888722258141- 10" 

10 0.8715278679388174731063544-10- 

1 1 0. 1979783472020383084286900- 10" 

12 0.4103178180353125619414324-10- 
Equilibration of the local extrema with the 

to 5.0- 10~ 18 : 

n b n 



1 .0000000000000000060373678 
0.999999999999997888979941 1 
0.5000000000001216148194572 
0.1666666666639271874501180 
0.4166666669859109153386033- 
0.8333333112815145481691497- 
0.1388889862738933258163839- 
0.1984098287973665146421103- 
0.2480734627092463176804164- 
0.2747848541489261879291146- 
0.2827881515524984459349078- 
0.2086709669366350082217004- 
0.3441995330913567239602395- 
following coefficients reduces this error 



10" 


1 


10" 


-2 


10" 


2 


10" 


-3 


10" 


-4 


10" 


-5 


10" 


-6 


10" 


-7 


10" 


-8 



1.7533876543770903961757996 

1 0.8503916537808109674449984 

2 0.1052086936309369262175803 

3 0.8722104733315565035195027-10- 

4 0.5434368311501568988037582-10- 

5 0.2711543491306964588901369-10" 

6 0.1128132888783054546376918-10- 

7 0.4024558229970401854905218-10- 

8 0.1256584419307581158766302-10" 

9 0.3488091466142981067488685-10- 

10 0.8715288225355426665019433-10- 

11 0.1979820112783685973416909-10" 

12 0.4092071997914099904014169-10" 



1.0000000000000000049913878 
0.9999999999999982006556032 
0.5000000000001063630793784 
0.1666666666642173677701902 
0.4166666669575619866774579 
0.8333333129071360606767972 
0.1388889803905621871712292 
0.1984099684263107292542263 
0.2480712594971168345247335 
0.2748077480706561519632930 
0.2826376902452534679601836 
0.2092376435267840500466468 
0.3432678789827820176761249 



10" 


1 


10" 


2 


10" 


2 


10" 


3 


10" 


4 


10" 


5 


10" 


6 


10" 


7 


10" 


8 



6. Summary 

Besides some generic algorithms to compute the Chebyshev series of inverse 
polynomials, there are two specific aspects that facilitate this task: (i) the expansion 
coefficients can be derived from the partial fractions of the inverse polynomial, (ii) 
Expansion coefficients with indices larger than the polynomial degree are recursively 
linked to those of lower order, (iii) An algorithm has been presented which derives 
a polynomial of a given degree such that the first terms of the Chebyshev expansion 
of the relative error of a given function represented by this polynomial vanish. 



Appendix A. Chebyshev Series of ln(l + x) 
The integral representation 

(A.l) H1 + x) -f£- 
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and term-by-term integration of (|2.1Q(I on the r.h.s. with (|2.9(l yield that the Cheby- 
shev coefficients of 

oo 

(A.2) f(x) = ln(l + x) = J2' UT*{x) < x < 1, 

n=0 

obey 

(A. 3) 2n/„ = a„+i,i(-3) - a„_i,i(-3), n > 1, 

explicitly [H| p. 88] QI] 

(A.4) /„ = K ' , n > 1 

from (|2.11|) . as tabulated in ^ 4.1.45]. The missing /o is 

, Atrt ,2 /- 1 ln(l+f) 3 + 2^ 

(A.5) fo = - , dx = 2\n , 

7T J y/x(l - x) 4 

because insertion of x = 1 in (|A.2(1 yields 

(A.6) /„ = 2 (/(!)-£/») 

and /(l) = In 2 and £~ =1 /„ = 21n(l + via 4.1.24]. 

Appendix B. Chebyshev Series of arctana; 
Integrating (|2.7|l over x with 

(B.l) / -dx ~ arctana: 

J l + x 2 

and US) we get p. 89] [Hj 

(B.2) arctana; = 2 > — ^ — ^T,(x), -1< x < 1, 

in particular at a; = 1 

(_)li/2j 



(B - 3 > I- S. 



From (IB. 2D and 11. 1511. the coefficients of 



=i^ 7 ,...^( 1 + ^2) J 



f-r, as arctanx 

(B.4) = 2^ g n T n {x) 



X n=0,2,4,6, 



as listed in (4.4.50)] 5 follow recursively, where go = 21n(l + V2) is obtained 
via 4.531.12]. 
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Appendix C. Chebyshev Series of arcsin: 



The series of arcsin x = X^ n=1 3 5 gnT n {x) starts with g\ = 4/ V. A combination 
of P (4.4.58)], (H3|), O and 2S (3.1. Id) 



(C.l) 

yields 
(C.2) 



4 °° 1 

- Y] -, — 2 T «( x ) 

7r l—vr 



(n + l)g n +i + (n- l)g r , 



8 n 



n = 2,4,6, 



7r n 2 — 1 

in this case, which can be unwound as g n = 4/ (7m 2 ). To find a formulation with 
controlled relative error, we would switch to h(x) — (arcsin x) jx = n —Q h n T n [x) 
to remove the zero in the spirit of example 15.51 With (|1.15|) , the expansion coeffi- 
cients are 



(C.3) 
(C.4) 



ho 

h-n+2 



= 8/?(2)/7T, 

= -K. 



7r(n + l) 2 ' 

where 0(2) 0.915965594177219015054603515 is Catalan's constant X Tab 23.3] O 
0.234.3]. 

Appendix D. Chebyshev Series of arcsin(x/^/2) 
The coefficients of 



(D.l) 



arcsin(x/v2) = 



oo 

£ 

n=l,3,5,... 



fori^n \Xj 



are found by partial integration of / arcsin([cos0]/\/2) cos(n9)d9 



(D.2) 



mr 



- cos[(rc-l)fl] ^ _ r cos[(n+l)6] M 
V2 - cos 2 9 J V2 - cos 2 9 



>> odd 



where 
(D.3) 
* cos(2m6>) 



d6 = G 2 



V2F(-±) « 2.6220575542921198104648395899, 



y/2[3F(^=) - 4J5(-4)] « 0.22577708482093539558499460534, m = 1, 



(D.4) 



G. s = 



o V2 - cos 2 u ^ v- L ^v^; 

are Complete Elliptic Integrals. To find a recurrence for these 

1 T s (x) 
i v /(2-x 2 )(l-x 2 ) ' 

we apply the method of 1, (17.1.4)] to the quartic y 2 = (2 — x 2 )(l — a; 2 ), with 
d{yT s {x))/dx = y(dT s (x)/dx)+T s ±(T 3 (x)-3T 1 (x)), insert O for the derivative 
on the r.h.s, replace the first y on the r.h.s. by y 2 /y = (T4/8— T 2 + 7/8)/y, expand all 
products with l|1.4|) . and finally insert the upper limit x — 1 where y{x)T s {x) = 0: 
(D.5) 



G s+ 3+G| s _ 3 | -3(G s+ i+G| s _!|) 



s-l 

■5 ? 

1 = 



[Gi 



8(G 



1+2 



Gi/_ 



|Z-2|, 



14G, 
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Inserting s = 1, 3, 5 and 7, for example, yields 



(D.6) 
(D.7) 
(D.8) 
(D.9) 



and generates fei to kg in l|D.2Jl from Go and Gi shown in l|D.3jl . A slowly converging 
series expansion is also known 0j 806.01]. With (|1.15(1 we find the coefficients 
/„ = 2fc„_i - f n -2 for 



3G 4 — 12G 2 ~\~ Gq 

5G 6 - 27G 4 + 15G 2 - G 
7G 8 - 41G 6 + 29G 4 - 3G 2 
9Gi - 55G 8 + 43G 6 - 5G 4 



0, 
0, 
0, 
0, 



(D.10) 

starting at 
fo = ■ 

(D.ll): 



1 x 
— arcsm —= 
x V2 



^ ' fn,T n (x), 



■1 < x < 1, 



arcsm 



V2 



x 



VT 



dx 



2 7T ^ 4q - 1 L 
9=1 y 



V2E 

;=o 

4q-2 ~ G± q 



2'(2Z + 1)[(2Z)!!] 2 

1.4866664932871034689603296833. 



The four coefficients a t that span 
(D.12) 



"-n— 1 — 

7T 



J n — 

7T 



+ (-) [ " /2l /o, 



start as follows: 



n 


a.\ 


«2 


a 3 


«4 


2 


-2 


4 


-4 


8 


4 


-26/9 


4 


-16/9 





6 


-638/75 


292/25 


-3428/225 


584/25 


8 


-22702/735 


212/5 


-513088/11025 


1536/25 


10 


-23722/189 


4652/27 


-6763436/33075 


191128/675 


12 


-463174/847 


2252/3 


-3558618544/4002075 


822272/675 


14 


-162508858/65065 


8691484/2535 


-2777152623884/676350675 


643269592/114075 



Because T 2 j(x) 
to a factor y/2: 



T*(x 2 ), the following numbers coincide with P (4.4.51)] up 



n 


fn 


n 


fn 


2 


0.3885303371652290716432228-10" 


-i 








4 


0.2885441422084471126676825-10" 


-2 


6 


0.2884218334475536563483289-10" 


-3 


8 


0.3322367192785279209254231-10" 


-4 


10 


0.4158477878052832866177270-10" 


-5 


12 


0.5496504525974164467345493-10" 


-6 


14 


0.7550078449371525934251585-10" 


-7 


16 


0. 1067193805629843129424091-10" 


-7 


18 


0.1542180379281470021561106-10" 


-8 


20 


0.2268114598545151963877153-10" 


-9 


22 


0.3383885639342775871004709-10- 


10 


24 


0.5108937524377197224216916-10^ 


11 


26 


0.7791139213632464421446539-10- 


12 


28 


0. 1198378589352895337866326- 10- 


12 


30 


0. 1856972621821342234640637- 10" 


13 


32 


0.2896189154386304361020997-lQ- 


14 


34 


0.4542792886328823081478511-10- 


15 


36 


0.7161678029265506176831289-10- 


16 


38 


0.1134144256904559996509711-10- 


16 
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Appendix E. Chebyshev Series of tp(x + 2) 
An expansion of the Digamma function PI] is P3 (6.3.16)] 

(E.l) ^ 2 + x ) = l-j + x Y /Tr -_ 



k=2 



k{x + k) ' 



where 7 w 0.5772 is Euler's constant. Employing a n ,i(— k) of l|2.4|l . 

(E.2) - = = - V \ ; TJx), -Kx<l, 



(E.3) 0(2 + ,) _ 1 - 7 + 2xg £ _^=_T„ W . 

The auxiliary definition 

00 

(E.4) K n = V - n = 0,l,2,... 

turns 1E.3|) with the aid of l|1.4|l into 

00 

(E.5) i>{x + 2) = (l- 7 -A'i)To(a:)-5^(-) n (lf n _i+if n+1 )r n (a;) ) -1 < x < 1, 

n=l 

where 

OO ^ 

(E.6) i^ n + i4T„ +2 =2V -== ; , n = 0,1,2,... 



Alternatives to the slowly converging original series l|E.4|) at small n are obtained 
in terms of the Riemann Zeta function £ after reducing the fraction in (|E.4|I and/or 



(IR6I) by fc - Vfc 2 - 1, 

M = Ei(i-^" /2 =E(-)^"y 2 )[c(2/ + 2)- 1 ], 



fc 2 V k 2 J ^ ' V I 

k=2 v 7 z=o v 



fe=2 ,„ v ^ J=1 

OO 



(E.9) X„ + if 2 = 2£(-)'( 1/2 )[C(20-1], 

1=1 ^ ' 



OO 



( E .io)A 1+ if 3 = 2j2(-y(( ) /2 ) + ( 1/ , 2 ))m-i)-i], 



00 n+l 

(E.l% + if„ +2 = 2 £ (_)i[ C ( 2 l-n)-l]£ 

/=[(n+3)/2] s=0 



n+fj /(a-l)/2> 



22 



RICHARD J. MATHAR 



n 


K n 


n 


K n 





0.6942240199692270653811973 




1 


0.1181923495113155830503315 




2 


0.2615575442260127035429158 


10" 


- 1 


3 


f\ /"* C\ r~ T A A i"\ A A — T A A AAA 1 A yi ( A r~ >—* / \ A A 

0.6357242927298094244957032 


io- 


- 2 


4 


0.1613702909326556648518537 


10" 


- 2 


5 


A J 1 i ai /. /". a j -J r* -i A A A T A / i it A A r~ 

0.4189942166841513997803225 


io- 


- 3 


t> 



0.1101726048982138724638504 


10" 


-3 


7 


0.2918277395837793537278094 


io- 


- 4 


o 

8 


0.7763995103341854698876680 


io- 


-5 


9 


0.2071 120322602199079344235 


io- 


-5 


10 


0.5534045978754736410165904 


io- 


-6 


11 


0. 1480224417758054637706871 


io- 


-6 


12 


0.3961806941781982189370558 


io- 


-7 


13 


0.1060807013890109056491206 


io- 


-7 


14 


0.2841134565373781348071928 


io- 


-8 


15 


0.7610594780500236739360477 


io- 


-9 


16 


0.2038876075855356359426642 


io- 


-9 


17 


0.5462507275750785409310247- 


io- 


10 


18 


0.1463563991421891531670298- 


io- 


10 


19 


0.3921418684181661587649434- 


io- 


11 


20 


0.1050708536652289553110610- 


io- 


11 


21 


0.2815309431326615497301838- 


io- 


12 


22 


0.7543503527420115661214215- 


io- 


13 


23 


0.2021259323594124792356794- 


io- 


13 


24 


0.5415919982188370035264436- 


io- 


14 


25 


0. 1451 186573479985220441655- 


io- 


14 


26 


0.3888434449455691142432431- 


io- 


15 


27 


0.1041901454404881356692972- 


10- 


15 


28 


0.2791764103483896329393674- 


io- 


16 











As a by-product, insertion of x = ±1 in (|E.3|) with ^(l) = — 7 shows 

k - 1 + Vfc 2 - 1 



(E.12) 



(E.13) 



E 

oo 

E 

k=2 



1 



kVk 2 
1 



l k + 1 + Vfc 2 - l 

fc + 1 + \/k 2 - 1 



ky/k 2 - 1 fc - 1 + 



1 



Linear combinations of these two equations are 

(K14) x (k~m+vw—i) 



(E.15) 



E 



i[(fc + VI 2 ^) 2 - 1] 
k + Vk 2 -l 



1 

2' 

= 1. 



= 0, 



, =a fcVF^Tp + VF^T) 2 - l] 

and these two can be combined to 



(E.16) 



OO 

E 



fc = 2 V^TKfc + Vk^l) 2 - 1] 
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